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Abstract 

The present parametric study deals with the meso-scale modelling of concrete subjected to cyclic com- 
pression, which exhibits hysteresis loops during unloading and reloading. Concrete is idealised as a two- 
dimensional three-phase composite made of aggregates, mortar and interfacial transition zones (ITZs). 
The meso-scale modelling approach relies on the hypothesis that the hysteresis loops are caused by 
localised permanent displacements, which result in nonlinear fracture processes during unloading and 
reloading. A parametric study is carried out to investigate how aggregate density and size, amount of 
permanent displacements in the ITZ and the mortar, and the ITZ strength influence the hysteresis loops 
obtained with the meso-scale modelling approach. 



1 Introduction 

Concrete subjected to cyclic compressive loading is characterised by a strongly nonlinear response. Un- 
loading and reloading is accompanied by nonlinear fracture processes on the meso-scale, such as debonding 
and slip between inclusion and matrix, which result in hysteresis loops [TJ [21 El H] - Furthermore, repeated 
loading of concrete subjected to cyclic loading may lead to failure although the stress applied is below 
the strength obtained for monotonic loading. 

Many phenomcnological constitutive models were developed over the last 30 years, which can describe 
the hysteresis loops and the fatigue strength observed in experiments. Examples are constitutive models 
based on bounding surface models O [7] or special combinations of damage and plasticity [Sj . These 
phenomcnological constitutive models have only limited predictive capabilities since they are based on 
curve fitting and not on the mechanics of the underlying physical processes. Therefore, they cannot 
predict reliably the response of concrete outside the range of experimental results available. 

Recently, the present authors have proposed a damage-plasticity interface model to describe the response 
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of concrete subjected to cyclic loading [5]. In this approach concrete was idealised as a three phase- 
composite consisting of aggregates, mortar and interfacial transition zones (ITZs). A special lattice 
model was used to descretise the different phases [TUl E] . The meso-scale modelling of hysteresis loops 
and fatigue strength was based on the assumption that the nonlinear fracture processes during unloading 
and reloading are caused by localised permanent displacements. The amount of permanent displacements 
was controlled by a damage-plasticity interface model which was applied to the mortar and ITZ phase. 
Aggregates were assumed to respond elastically. This meso-scale modelling approach may have better 
predictive qualities than macroscopic phenomenological models, since the macroscopic response can be 
designed by changing the response of the meso-scale constituents. Nevertheless, each of the three phases 
has a complex heterogeneous microstructure, which again requires numerical modelling by means of 
phenomenological models or multiscale analysis in which the mechanical response of the three phases is 
determined by their micro-structures [TJl E3 03] • 

In the present work, this meso-scale modelling direction is further pursued and a parametric study is 
performed to investigate the influence of several model parameters on the description of hysteresis loops. 
The first two parameters are the volume fraction of aggregates and the aggregate size, which are related 
to the random discretisation of the composite. Furthermore, three material parameters, which influence 
the amount of permanent displacements in the mortar and the ITZ and control the ratio of the strength 
of mortar and interfacial transition zone, are investigated. This parametric study is similar to other 
parametric studies presented in the literature on the modelling of concrete |15[ 116] . Nevertheless, the 
present work is focused on the meso-scale modelling of concrete subjected to cyclic loading. This is a 
direction which has not been pursued before to the authors knowledge. 



2 Modelling approach 

In the present study the response of concrete subjected to cyclic loading is analysed with a damage- 
plasticity interface model, which was developed in 0. In the following paragraphs, the modelling approach 
is briefly reviewed to introduce the parameters, which are investigated in the present study. The interface 
model relies on a combination of a plasticity model formulated in the effective stress space and isotropic 
damage mechanics. For the two-dimensional version of the model, the displacement jump u c = (it n , u s ) T 
at an interface is considered, which is transferred into strains e = (e n ,e s ) by means of the length h as 



The strains are related to the nominal stress <x = (er n , er s ) as 



ui is the damage variable, D is the elastic stiffness, e p = (e pn , e ps ) is the plastic strain and er = (cr n , <7 S ) 
is the effective stress. The elastic stiffness is 



where E and 7 are model parameters controlling both the Young's modulus and Poisson's ratio of the 
material. 

The small strain plasticity model is based on the effective stress <r and consists of the yield surface, 
flow rule, evolution law for the hardening parameter and loading and unloading conditions. A detailed 
description of the components of the plasticity model is presented in The initial yield surface is 
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Figure 1: Stress-strain response for fluctuating normal strains for fi — 1 (solid line) and fi — (dashed 
line). 



determined by the tensile strength f t , by the ratio s of the shear and tensile strength, and the ratio c of 
the compressive and tensile strength. The evolution of the yield surface during hardening is controlled by 
the model parameter fj,, which is defined as the ratio of permanent and reversible inelastic displacements. 
The scalar damage part is chosen so that linear stress inelastic displacement laws for pure tension and 
compression are obtained, which are characterised by the fracture energies Gft and Gf c . 

The eight model parameters E, 7, / t , s, c, Gft, Gf c and fi can be determined from a tensile, shear and 
compressive test of the material phase. 

The constitutive response of the interface model is demonstrated by the stress-strain response for fluctu- 
ating normal strains for fi = 1 and fi — (Figure [1]). The normal strain is increased to point A (A ). Then 
the strain is reduced to point B (B ) and again increased to point G (G ). The parameter /1 controls the 
amount of plastic strains. For fj, = a pure damage- mechanics response is obtained and the stress-strain 
curve is unloaded to the origin. For fi = 1, on the other hand, a pure plasticity model is obtained. 
The unloading is elastic and the compressive strength is reached sooner than for /i = 0. However, the 
magnitude of the compressive strength is similar for \x — 1 and [i = 0, since the plasticity model is based 
on the effective stress. The constitutive model describes the loss of stiffness and permanent strains, but 
it does not represent the hysteresis loops for unloading and reloading. Instead, these loops are obtained 
from the structural analysis, which is discussed in the following section. 

The interface model is applied to the two-dimensional plane stress meso-scale analysis of concrete. A 
special lattice-type model developed by [TTj [TSJ [TU] is used to discretise the domain. The domain is 
decomposed into polygons by means of the Voronoi tessellation [H5] (Figured)- The lattice elements 
connect the centroids of the polygons. Each node possesses two translations and one rotation shown in 
the local coordinate system in Figure (2b- The degrees of freedom u c = {u±, v\, <fii, 112, V2, 02 } T of the 
two nodes of the lattice element are related to the displacement discontinuities u c = {u Cl v c } at the 
mid-point G of the interface by 

u c = Bu c (4) 



where 



B = 



-10 e 1 -e 
3 -1 -h/2 1 -h/2 



(5) 



In Eq. ((5|), the variable h is the element length and e is the eccentricity, defined as the distance between 
the mid-points of the lattice element and the corresponding polygon facet, respectively (Figure ^ . The 
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displacements u c at the point C are transformed into strains s — {e n ,s s } — u c /h (Figure [2b) . The 
strain e is related to the stress <r = {a n , a s } T by means of the two-dimensional version (at — 0) of the 
damage-plasticity interface model presented earlier. The element stiffness is 



K c = 4b t DB 

h 



(G) 



where A is the length of the facet. 



The nodes of the lattice elements are placed sequentially in the domain to be analysed. The coordinates of 
each node are determined randomly and a minimum distance d nl = 3 mm is enforced iteratively between 
the nodes [20 . For this iterative process, the number of vertices n for a chosen domain ^4d and the 
minimum distance d m determine the distribution of lengths of lattice elements. This relationship can be 
expressed in the form of a density 

nd 2 

P=^ (7) 

The resulting irregular arrangement of lattice elements is more suitable for fracture analysis, since the 
fracture patterns obtained are less sensitive to the arrangement of the lattice elements. 

The meso-structure of the heterogeneous material is discretised by lattice elements perpendicular to 
the boundary between the cylindrical inclusions and the matrix. The diameters of the inclusions are 
determined randomly from the cumulative distribution function 



Pf 
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(8) 



where Pf is the probability of the occurrence of the inclusion diameter d. Furthermore, a = </> m ax/0min, 
where m ; n and </> max are the minimum and maximum diameters, respectively [21] . A pseudo-random 
number generator is used to generate probabilities from which the diameter d is determined by means 
of the inverse of Eq. l[5]). This procedure is repeated until the chosen volume fraction p a of circular 
inclusions is reached. The inclusions are placed sequentially by means of randomly generated coordinates 
within the area of the specimen. For each set of generated coordinates, it is checked that no overlap with 
previously placed inclusions occurs. However, overlap with boundaries is permitted, i.e. no boundary 
effect is considered. 




i- 



c. 



V-2 i 



Ul 



02. 



U-> 



(6) 



Figure 2: Discretisation: (a) Lattice elements (dashed lines) connecting centroids obtained with the 
Voronoi tesselation (solid lines), (b) Degrees of freedom u\, Vi, <j>i, «2, «2 and (f>2 of the lattice element 
of length h in the local coordinate system. The point C at which the interface model is evaluated is in 
the center of the polygon facet at a distance e from the center of the lattice element. 
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Figure 3: Concrete in uniaxial compression: a) Geometry and loading setup, (b) Lattice (c) Voronoi 
polygons (d) Boundaries of the inclusions. 



Table 1: Model parameters. 

Phase E [GPa] 7 f t [MPa] s c G {t [J/mm 2 ] G fc [J/mm 2 ] p 

Mortar 30 0.33 5 2 20 400 40000 1 

Interface 48 0.33 1 2 20 80 8000 1 

Aggregate 120 0.33 - - - - 



3 Parametric study 

The present parametric study deals with the modelling of concrete subjected to cyclic loading. A three- 
dimensional concrete cylinder subjected to cyclic uniaxial compression tested by Sinha, Gerstle and Tulin 
[1] was idealised by means of a two-dimensional meso-scale plane stress model. Concrete was represented 
by a three-phase composite consisting of mortar, aggregates and interfacial transition zones between the 
two phases. 

The meso-structure was discretised by means of a random lattice based on the Voronoi tesselation of the 
domain. The mesh was generated with a vertex density of p = 0.6 and a minimum distance of d m — 3 mm. 
This vertex density is close to the saturated vertex density which was determined by Bolander and Saito 
in [lOj as p = 0.68. The aggregate distribution was obtained with <i max = 32 mm, <i m i n = 10 mm and 
p a = 0.3. The geometry, loading setup and one exemplary mesh are shown in Figure [31 

The concrete specimen subjected to monotonic and cyclic compression was analysed with the model 
parameters shown in Table [T] The axial average stress was determined as a — P '/ d, where P is the axial 
force, Au is the axial displacement and d is the width of the specimen. Correspondingly, the axial average 
strain was defined as e = Au/L, where L is the height of the specimen. The analysis results in the form of 
the average stress-strain curves for monotonic and cyclic loading are compared to the experimental results 
in Figure [U A reasonable qualitative agreement between the modelling response and the experimental 
results is obtained. The deformation pattern in the post-peak regime are presented in Figure [51 
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Figure 4: Stress-strain responses obtained from (a) the meso-scale analysis and from (b) experiments 
reported in pQ. 
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Figure 5: Deformations magnified by a factor of five shown for the monotonic loading. 



Five parameters, which are likely to have a strong influence on the modelling of concrete subjected to 
cyclic loading, are studied: 

1. Volume fraction of aggregates. 

2. Aggregate size. 

3. Ratio of permanent and total inelastic displacements in the ITZ phase. 

4. Ratio of permanent and total inelastic displacements in the mortar phase. 

5. Ratio of the mortar and interface strength. 

The first two parameters, which are related to the aggregate distribution, require random lattice gener- 
ations. Random meshes and the random positions of the aggregates influence the average stress-strain 
response. Therefore, the influence of different discretisations on the response of monotonic and cyclic 
compressive loading is initially analysed by means of five analyses with the same aggregate size and 
distribution, but with random geometry. The five different meshes are presented in Figure El 

The response for both monotonic and cyclic compressive loading was analysed. For cyclic loading, the 
specimens were subjected to four cycles. The strain was increased in steps to e m = 0.5%, 1%, 1.5%, 2%, 
with the strain being reversed after each step so that <r m = 0. The same loading history was applied for all 
the analyses in the present parametric study. The results for the influence of the random positions of the 
aggregates and the random orientation of the lattice elements are presented in Figure [7] for monotonic and 
in Figure[5]for cyclic loading in the form of average stress-strain relations. The random discretisation has 
a small influence on the pre-peak regime, whereby the peak and post-peak regime is strongly influenced. 
Nevertheless, the scatter of the simulations is reasonable compared to the one obtained experimentally 
by van Mier in [3] . 

The first of the five parameters is the volume fraction p a . The response for three volume fractions 
p a = 0.3, 0.15 and were analysed. For p & — the resulting, stress-strain curves exhibit snap-back in the 
post-peak regime. To avoid numerical problems, the tensile strength of the mortar phase was reduced to 
ft = 2.1 MPa. The other model parameters are the same as in Table [TJ The distribution of the randomly 
placed aggregates is presented in Figure [5] 

The average stress-strain response for the different aggregate volume fractions is presented in Fig- 
ures [TO] and [Tl] for monotonic and cyclic loading, respectively. For the monotonic response, a decrease 
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Figure 6: Five random meshes with the same aggregate density and size distribution. Lattice elements 
which describe the ITZs between aggregates and mortar are not shown to improve the clarity of the 
illustration. 
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Figure 7: Stress-strain responses obtained for monotonic loading for five meshes with different aggregate 
arrangements and discretisations. 



8 




-30 



-3 -2 
strain e [mm/m] 



Figure 8: Stress-strain responses obtained for cyclic loading for five meshes with different aggregate 
arrangements and discretisations. 
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(a) (b) (c) 



Figure 9: Lattices for aggregate densities of (a) p a = 0.3, (b) p a = 0.15, (c) p a = 0. Lattice elements 
which describe the interface between aggregates and mortar are not shown to improve the clarity of the 
illustration. 
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Figure 10: Stress-strain responses for monotonic loading for three aggregate densities of p a = 0.3, 0.15 
and 0. The tensile strength of the mortar phase for the specimen with p a = was reduced to ft — 2.1 MPa 
to avoid snap-back in the post-peak regime. 
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Figure 11: Stress-strain responses for cyclic loading for three aggregates densities p a = 0.3, 0.15 and 0. 
The tensile strength of the mortar phase for the specimen with p a was reduced to ft — 2.1 MPa to avoid 
snap-back in the post-peak regime. 
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(a) (b) (c) 

Figure 12: Lattice elements for the size range of (a) d ma x/rfmm = 3.2, (b) dmax/^min = 2.1, (c) 

^max/^min 1.3. 

of the aggregate fraction results in an increase of the compressive strength. This strong influence of 
the aggregate fraction is explained by the weak interfaces. With an increasing number of interfaces, the 
strength of the composite is reduced. For the cyclic response, a decrease of the aggregate volume fraction 
leads to a decrease of the size of the hysteresis loops. The occurrence of hysteresis loops is based on the 
existence of weak zones, in which permanent displacements localise. If the number of these localised zones 
of permanent displacements is reduced, the nonlinear fracture processes during unloading and reloading 
is decreased and the resulting hysteresis loops are smaller. 

The second parameter studied is the ratio of the maximum and minimum aggregate size d max / d m i n . Three 
size ratios of <i m ax/dmin = 3.2, 2.1 and 1.3 were chosen. For all three cases, the maximum aggregate size 
was d max = 32 mm. The meshes are displayed in Figure [T2l 

The average stress-strain curves for monotonic and cyclic compressive loading for the different aggregate 
size ratios is shown in Figures fT3l and fT4l The difference between the curves is small and falls in the range 
of scatters obtained by the initial study of the influence of the discretisation. Thus, the size ratio does 
not influence the monotonic and cyclic response strongly. The remaining three parameters are related 
to the material response and are independent of the discretisation used. Therefore, the results for these 
parameters are obtained with the same discretisation. 

The next two parameters studied are related to the ratio of permanent and total inelastic displacements in 
the ITZ and mortar phase. This ratio is equal to the model parameter [i. The present meso-scale approach 
for cyclic loading is based on the hypothesis that the occurrence of localised permanent displacements 
is the main cause for the presence of the hysteresis loops. Consequently, a reduction of the permanent 
displacements (decrease of /i) should result in smaller hysteresis loops. The ratios fi — 1, 0.95 and 0.9 
were used. The monotonic and cyclic stress-strain curves for the three ratios of the ITZ phase are shown 
in Figure [TBI and Figure [T5I respectively. The size of the hysteresis loops decreases with increasing ratio 
of permanent and total inelastic displacements. This corresponds to the hypothesis that the hysteresis 
loops are caused by localised permanent displacement. 

Furthermore, the influence of the permanent displacements in the mortar phase is studied. Again, the 
same values of fi = 1, 0.95, 0.9 were used. The average stress-strain relations for monotonic and cyclic 
loading are presented in Figures [17] and [18] The ratio of the mortar phase has almost no influence on 
the pre-peak regime, since almost all nonlinearities result from the nonlinearities at the ITZs. However, 
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Figure 13: Stress-strain responses for monotonic loading for three aggregate-size ranges of d max /d n 
3.2, d maK /d m [ n — 2.1 and d mSL ^/d m i n = 1.3. 
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Figure 14: Stress-strain responses for cyclic loading for three aggregate-size ranges of d mlix /d n 
^max/ ^min = 2.1 and d maK /d m [ n — 1.3. 
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Figure 15: Stress-strain responses for monotonic loading for \i — 1, 0.95 and 0.9 of the ITZ phase. 




-5 -4 -3 -2 -1 



strain e [mm/m] 

Figure 16: Stress-strain responses for cyclic loading for fi = 1, 0.95 and 0.9 of the ITZ phase. 
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Figure 17: Stress-strain responses for monotonic loading for fi = 1, 0.95 and 0.9 of the mortar phase, 
in the post-peak regime nonlinear displacements occur in the mortar phase. 

Finally, the influence of the ratio of the strength of ITZ and mortar was investigated. Three strength 
ratios of /i// m = 0.1, /i// m = 0.2 and /i// m = 0.3 were considered. The average stress-strain relations 
for monotonic and cyclic loading are shown in Figures [H] and [30] 

The strength ratio of ITZ and mortar has a strong influence on the monotonic and cyclic response. The 
compressive strength increases with an increased ITZ/mortar strength ratio, since the strength of the 
weak zones (ITZs) is increased. For the cyclic response, the size of the hysteresis loops decreases with 
increasing interface strength, since less permanent displacements take place at the interfaces between 
mortar and aggregates. 



4 Conclusions 

In the present work a meso-scale modelling approach developed in [0] was used for a parametric study 
of the modelling of concrete subjected to cyclic compression. Concrete was idealised by a three-phase 
composite consisting of mortar, aggregates and interfacial transition zones. Five parameters were investi- 
gated: aggregate volume fraction, aggregate size, ratio of the permanent and total inelastic displacements 
in the ITZ and mortar phase ,and the ratio of the strength of the interfacial transition zone and the mortar 
phase. The influences of the different parameters are summarised as following: 

• A decrease of the volume fraction of aggregate increases the ultimate strength of concrete, since 
fewer aggregates correspond to a reduction of ITZs, which weaken the material. A decrease of 
volume fraction reduces also the size of the hysteresis loops, since the number of ITZs, at which 
permanent displacements take place, are reduced. Localised permanent displacements are the main 
reason for the occurrence of hysteresis loops in the present meso-scale modelling approach. 

• The size range of aggregates does not strongly influence the response of concrete subjected to 
monotonic and cyclic loading. 
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Figure 18: Stress-strain responses for cyclic loading for \i = 1, 0.95 and 0.9 of the mortar phase. 
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Figure 19: Stress-strain responses for monotonic loading for interface-mortar strength ratios of /ti//tm 
0.1, 0.2 and 0.3. 
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Figure 20: Stress-strain responses for cyclic loading for interface-mortar strength ratios of fa/ftm = 0.1, 
0.2 and 0.3. 



• A decrease of the amount of permanent displacements in the ITZ phase results in a reduction of 
the size of the hysteresis loops. 

• A decrease of the amount of permanent displacements in the mortar phase results in a decrease of 
the size of the loops in the post-peak regime, in which localised permanent displacements occur in 
the mortar phase. 

• Finally, the ratio of the strength of the interfacial transition zone and the mortar has a strong 
influence on both the monotonic and cyclic response. An increase of the strength of ITZ leads to 
an increase of the compressive strength, since the amount of permanent displacements in the ITZs 
is reduced. Furthermore, the size of the hysteresis loops is decreased as well. 
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